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IHTRODUCarrO!! 


A  atethod  has  been  dereloped  for  the  deternlnatlcp  of  a  ccx^letc;  aet 
of  orbital  paraawtera  froa  a  fev  minutes  of  Doppler  data  recorded  In  the 
course  of  a  single  pass  of  a  satellite.  The  source  of  the  signal  may  be 
a  transmitter  In  the  satellite  or  a  ground-based  tiansaltter  reflecting 
a  algoaJ.  frca  the  satellite.  ^Hie  latter  transmitting  system  reqvilres 
more  costly  and  complex  equliatent  but  offers  reliability,  an  accurately 
Imown  transmitter  frequency,  and  a  stronger  geometry  for  a  van  accurate 
orbital  computation  vhen  the  number  of  receiving  stations  Is  limited. 

Since  It  was  desired  to  develop  a  rapid,  reliable,  and  moderately 
accurate  method  of  determining  the  orbital  parameters  of  a  satellite 
tracked  by  a  Doppler  system  employing  a  minimum  of  receiving  stations, 
emphasis  was  placed  on  the  development  of  a  solution  from  single  pass 
observations  recorded  at  froa  one  to  three  receiving  sites.  The  single 
pass  limitation  was  consider^  to  present  a  challenging  and  worthy 
problem  for  vblch  there  would  be  numerous  applications  If  a  reasonable 
solutloo  coidd  be  developed. 

In  the  past,  Doppler  data  have  been  used  primarily  to  measure  the 
slope  and  time  of  Inflection  of  the  frequency- time  cxirve  to  obtain  slant 
range  and  tlme-of-closest-approacb  Information.  This  Is  considered  to 
be  only  an  elemental  use  of  the  Information  In  the  Doppler  data.  Single 
pass  observations  from  one  receiver  have  been  demonstrated  to  contain 
sufficient  Information  for  satellite  orbital  determinations  of  8uffl'''?»nt 
accuracy  for  many  applications. 

Tor  exai^le,  it  may  be  desirable  to  know  the  orbital  parameters  as 
quickly  as  possible  after  launching  a  satellite.  The  orbital  parameters 
of  a  newly  launched  satellite  could  be  computed  within  mluates  after  the 
beginning  of  Its  free  flight.  Again,  after  atteaptlng  to  deflect  or 
steer  a  satellite  Into  a  different  orbit,  it  nay  be  desirable  to  kncv 
the  new  orbital  parameters  within  a  natter  of  minutes. 

The  fcUovlng  sections  will  discuss  *he  praetlcsllty  of  orbit  de¬ 
terminations  from  Doppler-  data  alone  and  will  Indicate  limitations  as 
veil  as  the  obvious  advantages  for  xhls  conceptually  slsqile  technique. 
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DESCREPTICN  OF  TRACKHTO  EQUUMENT  AND  MTA 


Doppler  obsermtlona  ounulat  of  recordings  of  Doppler  frequency.,  &s 
a  function  of  time.  Here  the  Doppler  frequency  is  defined  to  be  the  fre¬ 
quency  obtained  by  heterod^'nlng  a  locally  generated  signal  agtlnnt  the 
signal  received  from  the  satellite  followed  by  a  correction  for  the  fre¬ 
quency  bias  introduced  as  a  result  of  the  difference  between  the  frequency 
of  the  local  oscillator  and  that  of  the  signal  source.  Ih  t&ls  rejurt^ 
the  Doppler  frequency  is  defined  to  be  negative  when  the  satellite  Is 
approaching  the  receiving  site  and  positive  when  It  Is  receding.  If  the 
Doppler  frequency,  as  defined,  la  plotted  as  a  function  of  tiae,  one 
obtains  a  curve  of  the  form  shewn  In  Figure  1,  usually  referred  to  as  an 
"8"  •curve.  The  aoynuietry  of  the  curve  Is  typical  for  a  tracJclng  system 
with  a  ground-based  transmitter  and  a  receiver  separated  by  an  appreciable 
distance.  Only  for  a  satellite  whose  orbital  plane  hluecta  the  base  line 
will  the  Doppler  data  produce  a  syaaetrlcal  "B"  curre  with  a  reflection 
system.  With  a  satellite -home  transmitter,  the  "S"  curve  la  very  nearly 
symmetrical,  being  modified  slightly  by  the  Barth's  rotation  and  the 
refractive  effect  of  the  Ionosphere.  If  continuous  ohsezvatlans  are  made 
and  sampled  at  freqiient  Intermls,  such  as  or^f  per  second.  Figure  1  (a) 
Illustrates  an  analog  plot  of  the  data  avail.. jle  for  computer  input. 

However,  with  a  ground  based  transmitter  It  may  be  necessary  to  limit  the 
nusiber  of  observations  In  order  to  minimize  equlpnent  cost  and  coitP^exlty. 

For  example.  It  Is  possible  to  use  only  three  antenna  beams  and  provide 
three  sections  of  the  "3”  curve  as  shown  In  Figum  1  (b).  Another  possi¬ 
bility  la  the  use  of  a  scanning  emteona  been  to  provide  discreet  ^servatloos 
at  regular  Intervals  as  shown  In  Figure  1  (c).  Such  data  could  he  obtained 
by  an  antenna  with  a  thin,  fan-shaped  beam  which  scanned  the  sJ^  repeti¬ 
tively. 

Hie  data  In  any  of  the  forma  s-siggested  above  may  be  used  readily  as 
input  for  the  computing  procedure.  Whenever  possible,  tlus  Input  con¬ 
sists  of  the  ■t’otal  J>')ppler  cycle  count  over  a  variable  time  interval  re"  . 
than  the  Doppler  frequency  Itself  (l.e.  the  area  -aiuv-'  the  curres  or  euren 
of  curves  presented  In  Figure  1  (a)  and  1  (b)). 
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(c) 

Fig.  1 — Doppler  frequency-time  curves. 
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In  order  to  handle  the  Doppler  data  rapidly  and  "ccurately,  the 
Doppler  frequency  la  automatically  collated  and  dlgltlred  at  the  receiving 
altea.  Figure  2  ahowa  a  alnpllfled  block  diagram  of  a  DOPLOC  receiving 
ayotem.  Autoewitlc,  real-time  counting  of  the  Dpppler  frequency  requires 
a  algnal  of  high  quality,  that  la  one  with  very  small  randen  errors  Intro¬ 
duced  by  nolae.  Doppler  data,  which  are  essentially  noise  free,  are  made 
poaalble  In  the  DOPLOC  system  by  use  of  a  very  narrow  bandwidth,  phase- 
locked,  tracking  filter  (ref.  l)  following  the  receiver.  Significant 
taprovementa  in  the  slgnal-to-nolse  ratio  of  noisy  received  signals  are 
realized  by  extreme  reduction  of  the  system  bandwidth  throu^  the  use  of 
the  filter.  Bandwldths  adjustable  from  1  to  100  cyxOes  per  second  are 
available.  The  filter  Is  capable  of  phase- locked  operatlaa  when  a  signal 
la  as  weak  ns  36  decibels  below  the  noise,  (l.e.  a  nolse-to-slgnal  power 
ratio  oi.  kOOO).  The  filtering  action  Is  obtained  by  use  of  a  frequency- 
controlled  oscillator  that  Is  correlated  or  phase- locked  to  the  Input 
signal.  Cie  basic  block  diagram  of  the  tracking  servo  loop  Is  shown  In 
Figure  3*  !Q:aektng  Is  acconpUahed  with  an  electroole  servo  system  designed 
to  force  the  frequency-controlled  oscillator  to  follow  the  variations  of 
frequency  and  phase  of  the  Input  signal.  Correlation  Is  maintained  with 
respect  to  Input  signal  phase,  frequency,  first  time  derivative  of  frequency, 
and  with  a  finite  hut  small  phase  error,  the  second  time  derivative  of 
frequency.  This  Is  done  by  a  cross-correlation  detector  consisting  of  the 
phase  detector  and  filter,  or  equalizer  network.  Dtader  dynamic  conditions, 
the  control  voltage  to  the  oscillator  Is  so  filtered  In  the  equalizer  net¬ 
work  that  tracking  faithfully  reproduces  tbs  rate  of  change  of  the  Input 
frequency.  An  Inherent  feature  at  this  design  Is  an  effective  aceelexatlon 
memory  which  provides  smooth  tracking  and  extrapolation  thzmq^  signal 
dropouts.  Experience  with  signal  reception  froe  Earth  satellites  has  borne 
out  the  necessity  for  this  mmory  feetvre,  since  the  recelviid  signal  aaipll- 
tude  may  vary  widely  and  rapidly.  ISie  filter  work*  throu^  signal  null 
periods  very  effectively  without  losing  Irw'k.  jn  addition,  this  meao 
provides  effective  tracking  of  the  desired  Doppler  slgsal  In  the  pr^eenue 
of  lnterferl:ig  clgnals  when  severed.  s«telll'*'«B  are  within  receiving  nuge 
at  the  same  time. 
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STATION  USING  BRL  TRACKING  FILTER 

FIGURE  2 
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The  slgnal-to-nolse  power  Imprcveinent  furnished  by  the  tracKlng 
filter  Is  equal  to  the  ratio  of  the  .lr?put  source  noise,  hsocdwldth  to 
the  filter  bandwidth.  The  Internal  noise  generated  by  the  filter  Is 
negligibly  snail  at  all  bandwldths.  The  relation  between  Input  and 
output  slgnal-to-nolse  la  shown  In  Figures  k  and  5.  In  a  typical  case 
wl..h  a  receiver  bandwidth  of  10  kc  and  a  filter  bandwidth  of  5  ci>a,  a 
signal  burled  27  db  down  In  the  noise  will  appear  at  the  filter  outpit 
with  a  6  db  slgnal-to-nolse  ratio.  An  experimental  Investigation 
(ref.  2)  has  been  made  of  the  relation  between  slgnal-to-nolse  ratio 
and  the  uncertainty  or  random  error  In  measuring  the  frequency  of  a 
Doppler  signal.  The  test  results  showing  R.M.S.  frequency  error  as  a 
function  of  slgnal-to-nolse  ratio  and  tracking  filter  bandwldths  are 
shown  In  Figure  6.  For  the  example  cited  previously,  a  signal  27  db 
down  In  the  noise  can  be  read  to  an  accuracy  of  0.13  cps.  An  Integration 
time  or  counting  time  Interval  of  one  second  was  used  for  these  measure¬ 
ments  . 

The  tracking  filter  can  be  equipped  with  a  signal  search  and  autena- 
tlc  lock-on  system.  Signals  30  db  down  In  the  noise  at  typical  Doppler 
frequencies,  frem  2  to  Ilf  kc,  can  be  detected  within  a  fraction  of  a 
second  and  the  filter  phase- locked  to  the  signal.  With  this  equlinent, 
signal  acquisition  and  lock-on  have  beccaw  routine  In  field  operations. 

The  DQPLOC  system  has  been  used  extensively  for  satellite  tracking. 

The  Inherent  high  sensitivity  of  the  receiving  system  to  signals  of  very 

•20 

low  energy  (2  x  10  watts,  -  197  dbw,  or  0.001  microvolts  across  ohms 
for  a  threshold  signal  at  1  cps  bandwidth)  has  permitted  the  use  of  con¬ 
ventional,  low  gain,  wide  coverage,  antennas  to  achieve  horizon  to  hori¬ 
zon  tracking  at  great  ranges.  It  has  been  found  to  be  practical  to  change 
bandwldths  over  the  selectable  range  of  1  to  100  cps  In  accordance  with 
the  Infoimatlon  content  of  the  signal  and  thus  achieve  maximum  slgnal-to- 
nolse  ratio.  Since  the  key  to  successful  determination  of  orbits  lies  In 
obtaining  data  with  small  values  of  random  anl  systematic  error,  t  high 
quality  dfttii  output  of  the  DOPLOC  system  has  been  an  Importan'*'  feature. 

An  orbital  solution,  developed  specif! oa.'Ty  for  this  system,  has  yielded 
relatively  accurate  results  with  a  surprisingly  small  number  of  DOPLOC 
tracking  cbservatlons . 
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WAVEFORMS  VS  NOISE  TO  SIGNAL  RATIO 
FOR  I  CPS  BANDWIDTH 


OUTPUT  SWNAL  TO  NOISE  RATIO  -Sk 


THE  ORSmi  SO-UJT  ICH 


The  method  of  solution  conslats  of  a  curve-fitting  procedure,  in 
which  a  compatible  sat  of  approximations  for  the  orbital  parameters, 
are  improved  by  successive  differential  corrections.  The  latter  are 
obtained  from  a  least-squares  treatment  of  an  over-detenalned  system 
of  equations  of  condition.  The  imposed  limitation  of  single  pass 
detection  permits  several  assimipticnB  vdilch  considerably  simplify  the 
computing  procedure.  Among  these  is  the  assuntptlon  that  the  Earth 
may  be  treated  dynamically  as  a  sphere  while  geometrlcully  regarding 
it  as  an  ellipsoid.  In  addition,  it  is  assumed  that  no  serious  loss 
in  accuracy  will  result  if  drag  is  neglected  as  a  dynamic  force.  With 
these  assumptions,  it  is  apparent  that  the  satellite  nay  be  regarded 
as  movi.jg  in  a  Keplerlan  orbit.  An  additional  simplification  in  the 
reduction  of  the  tracldng  data  is  warranted  if  the  frequency  of  the 
system  exceeds  100  megacycles;  for  it  then  becomes  feasible  to  neg¬ 
lect  both  the  atmospheric  and  ionospheric  refraction  of  the  transmitted 
signal. 

In  formulating  the  problem  mathenatlcally,  it  is  helpful  to  regard 
the  instrumentation  as  an  Interf ercneter .  In  this  sense,  the  total 
number  of  Doppler  cycles  observed  within  any  time  Interval  will  provide  a 
measure  of  the  change  in  slant  range  frcm  the  receiver  to  the  satellite 
if  the  transmitter  is  air-borne,  or  in  the  sum  of  the  slant  ranges  from 
both  the  transmitter  and  receiver  to  the  satellite  if  the  signal  origi¬ 
nates  on  the  ground  and  is  either  reflected  or  retransmitted  by  the 
satellite.  Assuming  the  latter  for  the  discussion  which  follows,  let 

g  ( t )  be  defined  as  the  change  in  the  sum  of  the  two  slant  ranges . 

J 

It  follows  frum  Figure  7  that 

gj  (t^)  -  (TS2  +  HSg)  -  (TSj^  +  RS^)  ,  (1) 

where  T  is  the  position  of  the  transmitting  sit'.,  R,  the  location  01  tne 

1/ 

^tb  receiver,  S  ’'ce  position  of  the  satellite  at  time  t^ ,  and  B  the 
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mK>8L£ll  OCOMETRY 

FIGURE  7 
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position  a!-  time  gj  cliange  in  the  sum  of  the  nlaat 

ranges  from  the  satellixe  to  the  transrltter  and  to  the  Jth  receiver 
in  the  time  interval  from  t^^  to  t^.  It  is  worth  noting  that,  if 
this  time  interval  is  equal  to  one  second  and  >,  is  the  wavelength  of 
the  transmitted  signal,  ^  ^  equivalent  to  the  Doppler 

frequency  for  the  Jth  receiver  at  the  time,  (t^^  +  0.5  sec.). 

The  mathematical  development  of  the  computing  procedure  has  been 
presented  in  reference  3  and  will  not  be  repeated  here.  Rather,  we 
will  confine  our  remarks  to  a  summary  of  the  more  Important  phases  of 
the  method.  Hie  solution  consists  of  Inpmvlng  a  set  of  position  and 
velocity  components  which  have  been  approximated  for  a  specific  time. 

Hie  latter  will  be  defined  as  t^  and  in  general,  will  be  within  the 
time  iuterval  over  which  observations  have  been  recorded.  The  com¬ 
puting  procedure  ‘is  outlined  in  Figure  8.  initial  approxlaationa  for 
position  and  velocity  uniquely  define  a  Keplerlan  orbit  which  may  be 
described  in  te’-ms  of  the  following  orbital  -’ters: 

a  s  semi -major  axis, 
e  s.  eccentricity, 
a  a  mean  anomaly  at  epoch, 

1  s.  inclination, 

n  s  right  ascension  of  the  ascending  node, 

<0  a  argument  of  perigee. 

After  these  parameters  have  been  determined,  the  position  of  the 
satellite,  and  then  gj  (t),  may  readily  be  cotqwted  as  a  function  of 
time.  Ccmparlng  the  computed  values  of  gj  (t)  with  the  observed  values 
of  the  same  quantity  and  acoumlng  more  than  six  observations,  a  set 
of  differential  corrections  for  the  Initial  approxlmatlonr-  o-  position 
and  velocity  may  then  be  obtained  from  a  standard  least-squares  treatment 
of  the  resulting  over-determined  system  of  equations.  Hie  corrections 
are  applied  to  the  initial  approximations  and  the  casputation  la  lte..ated 
until  con/ergencs  is  achieved. 
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DIAGRAM 
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FIGURE  8 


Dlls  ccaaputlng  procedure  esset-isjJLy  determines  oaly  those  seg¬ 
ments  of  the  orbit  confined  within  the  Intervals  of  observation.  By 
constraining  the  satellite  to  Keplerlan  motion,  the  parameters  a,  e,  a, 

1,  n,  and  (D  are  likewise  determined  In  the  course  of  the  ccmputatlon; 
and  these  serve  to  provide  an  estimate  of  cotton  over  the  entire  orbit. 

On  the  other  hand.  It  has  been  found  Impractical  to  fit  an  entire  ellipse 
to  the  observations  by  solving  fcr  the  orbital  parameters  directly. 
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miTIAL  ORBmL  APP.’’COaK»TK»IS 


Convergence  of  the  computation  rests  prunarlly  t^oa  the  adequacy 

I  I 

of  the  Initial  approximations  for  position  and  velocity.  It  has  been 
establlbhed  that,  for  a  system  consisting  ci'  a  single  receiver  and  an. 
earth-bound  transmitter  at  opposite  ends  of  a  400  mile  base  line,  con¬ 
vergence  is  assured  when  the  error  in  each  coordinate  of  the  Initial 
estimate  Is  not  In  excess  of  50  to  75  miles  and  the  velocity  cootponents 
are  correct  to  vlthln  l/2  to  1  mile  per  second.  When  tte  signal  source 
Is  carried  by  the  satellite  a  unique  solution  la  Imiposslble  vlth  obser¬ 
vations  from  a  single  receiver.  Eovever,  If  single  pass  measurements 
are  available  from  two  or  more  receivers,  with  either  a  ground-based  or 
6U1  alr-bome  transmitter,  the  system  geometry  la  greatly  strengthened. 
Convergence  ofn  then  be  expected  when  the  initial  approximations  are, 
within  150  to  200  miles  of  the  correct  value  In  each  coordinate  and 
1  to  2  miles  per  second  In  each  velocity  component.  larger  errors  may 
occasionally  be  tolerated,  but  the  figures  presented  are  Intended  to 
specify  limits  within  which  convexgence  may  be  reasonably  assured. 

Tberefoire,  it  has  been  necessary  to  develop  a  siqpartlng  cosgnita- 
tlon  to  provide  relatively  accxirate  Initial  approxlaatlans  to  position 
and  velocity  for  the  primary  costputation.  Several  successful  methods 
have  been  developed  for  this  phase  of  the  problem;  but  discussion  will 
be  confined  to  a  few  applications  of  a  differential  equation  derived. 

In  reference  5,  to  approximately  relate  the  motion  of  the  satellite  to 
the  tracking  observations.  If  the  transmitter  Is  earth-bound,  this 
equation  Is  of  the  form 


where  the  slant  ranges  from  the  transmitter  and  the  receiver  to  the 
satellite  are  respectively  and  pj.  g^  Is  the  second  time  derivative 
of  the  function  defined  by  equation  (l).  In  deriving  equation  (2), 

It  was  asBimie!:'  that; 
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1)  la  angti.1  nr  measurement,  the  f!\tKlllt8  Is  within  ten  degrees 
of  the  Instnmientatlon  site, 

2)  the  Earth  Is  not  rotating, 

,3)  the  satellite  moves  la  a  circular  orbit. 

With  these  assumptions,  A  may  he  shown  to  be  approximately  equal  to 
vV(GS)  and  hence,  constant  for  a  circular  orbit  since  R  Is  the 
Earth's  radlxn:,  v  Is  the  velocity  of  the  satellite,  and  G  Is  the  mean 
gravitational  constant. 


The  first  application  to  be  considered  will  be  for  a  systeta  la 
which  the  transmitter  Is  carried  by  the  satellite.  For  the  receiver 
In  such  a  system,  equation  (2)  reduces  simply  to 


(3) 


If  measurements  of  the  rate  of  change  of  the  Dpppler  frequency,  f are 

J 

made  for  two  different  times,  t  and  t^,  and  Doppler  frequencies,  f,, 

O  A  J 

are  observed  at  regular  Intervals  between  t  and  t. ,  we  note  that 


X  fj  (t^), 

X  fj 

X  f ,  (t^). 


Pj  (t^)  "  X  fj 

Pj  (^i)  -  pj  + X  r 


fj  (t)  dt. 


where  X  Is  the  wavelength  of  the  signal  and,  (t^)  and  (t^)  are  tho 
only  unlmowns.  Combining  equations  (If)  with  equation  (3)  yields 
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vtalcb  vltb  the  last  relation  of  equations  (4)  deteialnes  elaat  range 
as  a  fimctlon  of  tlae.  Uiese  results  lauy  h;  used  with  equation.  (3)  to 
establish  a  valtie  for  A  fron  vhlch  an  excellent  asproxlaatloci  of  the 
velocity  of  the  satellite  aay  be  obtained,  ilo  additional  Infoimtlon  can 
be  extracted  when  observations  are  limited  to  those  from  a  single  receiver. 
However,  If  measurements  from  three  or  more  receivers  overlap  in  time, 
a  set  of  approximations  to  the  position  and  velocity  cosoonents  any  be 
determined  by  a  stral^tforward  trianculatlon  procedure.  When  data  from 
only  two  receivers  are  available,  an  estimate  of  position  and  velocity 
aay  still  be  obtained  for  a  time  idileh  lies  within  the  Interval  of  obser¬ 
vation  of  both  receivers.  If  the  results  of  the  cosvutatlon  described  by 
eqmtlon  {5)  are  combined  with  the  asaimptlon  of  circular  motion.  Par  an 
epoch  time,  selected  so  that  the  satellite  Is  near  tha  zenith  of  the 
Instruments tlun  site,  we  may  safely  assume  that  the  vertical  coaqKxient 
of  velocity  Is  small  and  can  well  be  approximated  by  zero,  using  the 
results  of  the  coaqiutlng  procedure  described  above,  slant  ranges  for  the 
epoch  time  aay  be  coaqnited  for  each  receiver;  and  In  the  process,  Sn 
estimate  for  the  velocity  of  the  satellite  will  be  obtained.  Combining 
these  three  results  with  the  Doppler  frequency  meastuements  from  the  two 
receivers  for  epoch  time,  we  msy  readily  determine  the  renalnlng  velocity 
components  and  all  three  position  eoo'*dlnates.  m  this  develppaeat,  no 
account  has  been  taken  of  the  difference  in  frequency  between  the  trsna- 
mltter  In  the  satellite  and  the  reference  oscillator  on  the  ground.  If 
both  are  stable,  a  constant  frequency  error,  or  bias,  will  be  Introduced. 

In  general,  this  error  is  so  large  that  it  must  be  corrected  before  applying 
the  above  procedure.  Most  methods,  for  determining  tha  bias,  asstme 
ayMKtry  ahcmt  the  Inf  lection  point  and  use  this  charaeterlatlc  of  the  "S" 
cuive  to  deteimlne  the  Inflection  time  as  accurately  as  possible.  Since 
the  latter  is  also  the  time  of  elcsest  approach  of  tha  aatelllre  to  the 
receiver,  the  Doppler  frequency  should  be  sero.  therefore,  the  bias  la 
simply  the  observed  frequency  at  the  Inflection  time. 


The  second  application  considers  a  system  In  which  the  trans¬ 
mitter  Is  earth-bound  so  that  the  signal  travels  frcn  the  Saith's 
surface  to  the  satellite  and  back  to  one  or  more  receivers  on  the 
Garth’s  surface.  For  this  problesj  equation  (2)  applies.  Let  us 
define  a  right-hand  rectangular  coordinate  systqm  as  shown  In  Figure 
9  with  the  origin  at  the  transmitter  and  the  Z-axla  positive  In  the 
direction  of  the  vertical,  nie  y-axls  Is  foimed  by  the  Intersection 
of  the  tangent  plane  at  the  transmitter  with  the  plane  determined  by 
the  transmitter,  the  J^th  receiver,  and  the  Earth's  center.  The  re¬ 
ceiver  will  then  be  at  the  known  point  (0,  yj,  Zj).  If  the  variable 
point  (x,  y,  z)  Indicates  the  position  of  the  satellite,  the  slant 
ranges  from  the  transmitter  and  the  receiver  are  respectively 
given  by 


from  which  It  follows  that 


In  the  three-beam  mode  of  operation,  the  satellite  will  be  approximately 
In  the  yz-plane  at  t^,  which  Is  defined  as  the  time  halfway  between  the 
Initiation  and  termination  of  tracking  In  the  center  beam.  Let  the 
satellite’s  position  and  velocity  at  this  time  be  defined  as  (x^,  y^, 
and  (x^,  y^,  i^),  respectively.  Obviously,  may  be  approximated  by  zero 
sod  as  before,  may  also  be  set  equal  to  zero.  Equations  (7)  then  reduce 
to 
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2l 


Let  and  be  the  Dopplei  frequency  and  rate  of  change  of  frequency 
for  the  ^th  receiver  at  t^.  It  follows  that 


^^To  '^Jo^ 


(9) 


From  equation  (2),  we  conclude 


^^o  ^Jo 


(10) 


Expressing  equations  (9)  and  (10)  In  terms  of  the  position  coordinates 
and  velocity  ccnponents  of  the  satellite  at  time,  t^,  yields 


Let  us  assuae  a  specific  orbital  Inclination.  With  our  previous 
asBioaptlon  of  circular  motion,  y^  may  readily  be  confuted  as  a  f'^irtlon 
of  y^  and  z^.  'ihon  equations  (U)  and  (12)  will  likewise  provide  fj^ 
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and  f.  aa  ftmctlons  of  position  In  the  yz- plane.  Thus,  for  a  given 

Jo 

Inclination,  families  of  curves  may  be  computed  and  plotted  In  the 

yz-plane  for  both  f,  and  f .  .  Figure  10  presents  such  a  plot,  for  an 
Jo  Jo 

Inclination  of  80  ,  with  the  transmitter  and  receiver  separated  by 

miles  and  with  both  located  35°  off  the  equator.  To  attain  synanetry 
and  simplify  the  construction  of  such  charts,  Zj  was  assumed  to  be  zero, 
which  Is  a  reasonable  approximation  for  this  approach  to  the  problem. 

If  similar  charts  are  prepared  for  a  nuaiber  of  Inclinations,  satisfactory 
Initial  approximations  may  be  rather  quickly  and  easily  obtained  by  the 
following  operations : 

1)  Assume  an  Inclination.  This,  of  course.  Is  equivalent  to 
selecting  a  chart.  Accuracy  is  not  essential  at  this  stage 
since  the  estimate  may  be  la  error  by  15°  without  preventing 
convergence . 

2)  Enter  the  chart  with  the  obeerved  vaU.ues  of  f,  and  f .  to 

(JO 

determine  aa  appropriate  position  within  the  yz-planc. 

3 )  Approximate  the  velocity  ccmponents .  Ihese  should  be  conslsten  t 
with  the  aaaumptlon  of  circular  motion,  the  height  determined 

In  step  2),  and  the  assuaed  Inclination. 

4)  Determine  the  position  and  velocity  eompooents  In  the  coordi¬ 
nate  system  for  the  primary  solution  by  an  apaprpprlate 
coordinate  tranefonaation. 

In  addition  to  the  graphical  method,  a  digital  solution  has  been 

devleed  for  equations  (11)  and  (12).  As  in  the  previous  developsient, 

we  have  two  measuremente  available  and  desire  to  detezmlne  three  unknoms. 

In  this  approach,  one  unknown  is  dstarmlned  by  establishing  an  upper  bound 

and  asaumlng  a  value  which  Is  a  fixed  distance  from  this  bound.  The 

distance  has  been  selected  to  place  the  variable  between  Its  upper  and 

lower  bounds  In  a  position  which  Is  favorable  for  convergence  of  the 

primary  computation.  In  this  method,  we  chose  to  start  by  approximating 

z  .  It  may  be  observed  in' Figure  10  that,  for  larger  values  of  f.  ,  the 
o  Jo 

maximum  value  o'*  z^  occurs  above  either  the  transmitter  or  receiver  while. 
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for  smaller  values  jf  f .  ,  the  maximum  value  of  z  occurs  over  the  mid- 

JO  o 

point  of  the  base  line.  Ihc  first  step  1".  the  coniputatlon  Is  to  determine 

a  maximum  '/alue  for  z^.  To  this  end,  Is  eliminated  froa  eqwtlona  (U) 

and  (12)  to  yield  an  expression  which  varies  only  In  y^  and  z^.  A  appears 

In  this  expression,  but  It  la  also  a  function  of  these  variables,  llie 

resulting  equation  may  be  solved  by  numerical  methods  for  z^  with  y  «  0 

and  then,  solved  a  second  time  for  with  y^  ■  l/2  y^.  The  larger  of 

these  results  Is  to  be  used  as  a  value  for  (z  )„  which  is  defined  to  be 

the  maximum  possible  value  of  z^.  Assuming  the  altitudes  of  all  satellites 

to  be  In  excess  of  75  miles,  we  may  ccmclude  frcm  the  general  characteristics 

of  the  family  of  curves  for  fj^  In  figure  10,  that  the  satellite’s  altitude 

will  differ  from  (*„)u  “o  more  than  100  miles.  Since  an  error  of  50 

miles  may  be  tolerated  in  the  approximation  for  each  coordinate, 

Qz^)|^  -  5^  I*'  a  suitable  value  for  z^.  With  the  altitude  thus  deteimlned, 

we  may  solve  equations  (U)  and  (12)  for  y^  and  y^.  In  the  process  A,  and 

hence  the  velocity,  will  be  detexmlned.  With  z^  assmed  as  zero,  x^  may  be 

readily  evaluated  to  complete  the  Initial  approximations  which  consist  of 

the  position  (0,  y^,  z^)  and  the  velocity  (x^,  y^,  O).  It  Is  worth  noting 

that  there  Is  a  pair  of  solutions  for  y^  and  y^.  further,  the  method  does 

not  determine  the  sign  of  x^.  If,  in  addition,  we  accept  the  possibility 

of  negative  altitudes  for  the  mathematical  model,  ve  arrive  at  eight 

possible  set?  of  Initial  conditions  which  are  approximately  symmetrical 

with  respect  to  the  base  line  and  Its  vertical  bisector.  It  Is  an  Interesting 

fact  that  adl  eight,  when  used  as  input  for  the  primary  ccmputatlon,  lead 

to  convergent  solutions  which  exhibit  the  saiss  type  of  syametiy  as  the 

zpproxlmatlons  theauelves.  Of  course.  It  Is  trivial  to  eliminate  the  fbur 

false  solutions  which  place  the  orbit  underground.  IVrther,  two  additional 

solutions  nay  be  eliminated  by  noting  that  the  order  in  which  the  satellite 

passes  through  the  three  antenna  beams  determines  the  sign  of  x  In  the 

o 

two  remaining  possibilities,  Is  observed  to  have  opposite  signs.  Since 
the  y-axls  of  the  DOPLOC  system  has  been  oriented  frea  west  to  east,  the 
final  ambiguity  may  be  resolved  by  aasuming  an  eastward  component  of 
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Telocity  for  the  aatellite  >  certelnly  a  valid  aaeiaqptlca  to  date.  Ih 
any  event,  aJJ.  ambiguity  may  be  reBoved  rrcm  the  aolutlon  by  the  addition 
of  one  other  receiver.  Moreover,  this  would  significantly  improve  the 
geoeutry  of  the  systes  and  thereby  strengthen  the  solution. 

Hie  first  method  presented  in  this  section  is  Intended  for  use  with 
a  satellite  which  carries  its  own  transmitter.  (Riese  data  are  generally 
recorded  contlnuoiisly  as  in  Figure  la.  nie  other  two  methods  have  been 
developed  for  a  system  which  provides  observations  of  the  type  displayed 
in  Figure  lb  where  the  signal  source  is  on  the  Earth’s  surface,  ibe  plot 
shown  in  Figure  Ic  la  also  for  a  system  with  an  earth-bound  transmitter; 
and  the  last  two  methods  may  be  applied  to  such  data  if  minor  modifications 
are  made  in  the  procedures.  Indeed,  with  any  tracking  systaa  that  provides 
observatloas  of  satellite  velocity  components,  eqimtlon  (2)  furnishes  an 
adeqijate  base  for  establishing  an  approoclmate  orbit  to  serve  as  an  initial 
solution  which  may  be  refined  by  more  sophisticated  methods. 
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RESUiars  0?  orbital  cckputaticm 

Niunerous  convergent  solutions  have  been  obtained  with  actual  field 
data  froB  a  system  consisting  of  a  transmitter  at  Fort  Sill,  Oklahoma, 
and  a  single  receiver  at  Forrest  City,  Arkansas.  Hils  system  eomiilex 
provides  a  base  line  of  1^3h  miles.  In  addition,  several  orbits  have  been 
established  from  field  data  for  satellites  which  carried  the  signal  source. 
For  the  latter  mode  of  operation,  receivers  were  available  at  both  Forrest 
City,  Arkansas,  and  Aberdeen  Proving  Ground,  Maryland,  providing  a  base 
Hn,.  of  863  miles.  Results  will  be  presented  for  both  types  of  tracking. 

Hie  Initial  successful  reduction  for  the  Fort  Sill,  Forrest  City 
system  vau  achieved  for  Revolution  9937  of  Sputnik  XU.  Hie  DOPLOC  obser¬ 
vations,  as  well  ae  the  results,  ace  presented  In  Figure  11.  Msasurements 
were  recorded  for  28  seconds  In  the  south  antenna  beam,  7  seeoods  in  the 
center  beam,  and  12  seconds  In  the  north  beam,  with  two  gaps  In  the  data 
of  75  seconds,  Hius,  observations  were  recorded  for  a  total  of  k7  seconds 
within  a  time  Interval  of  3  minutes  and  17  seconds.  Using  the  graphical 
method  described  In  the  previous  section  to  obtain  Initial  apprcadmatlons, 
convergence  was  achieved.  In  three  Iterations  on  the  first  pass  through 
the  computing  machine.  It  will  be  noted,  in  the  ccaparlscm  of  DQPLOC  and 
Space  Track  results,  that  there  Is  good  agreement  In  a,  e,  1,  and  n,  par¬ 
ticularly  for  the  latter  two.  HjIs  Is  characteristic  of  the  single  pass 
solution  when  the  eccentricity  is  small  and  the  ccqputatlonal  Input  Is 
limited  to  Doppler  frequency.  Since  the  orbit  Is  very  close  to  being 
circular,  both  a  and  a  are  difficult  for  either  the  DOPLOC  System  or  SPi^e 
Track  to  determine  accurately.  However,  ae  a  ree\d.t  of  the  aaall  eccen¬ 
tricity,  ((D  o)  Is  a  good  approximation  of  the  angular  distance  along  the 
orbit  from  the  equator  to  the  position  of  the  satellite  et  epoch  tlr»  ‘uid 
as  such  provides  a  basis  of  comparison  between  the  two  systMM.  A  ccaqparl- 
Bon  of  this  quantity  Is  Included  In  Figure  11.  To  summarize,  when  limited 
to  single  pass,  single-receiver  observations,  the  DCFLUC  System  provides 
an  excellent  deteivl nation  ■,T  the  orientation  of  the  orbital  plane,  a 
good  determination  of  the  shape  of  the  orbit,  and  e  ••  *.lr-to-poor  determi¬ 
nation  of  the  orientation  of  the  ellipse  within  the  orbital  plana. 
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Although  a  and  m  have  been  accurately  beteralned  cn  occasion,  the 
InterljR  DOTLOC  system  with  Its  limitations  falls  to  provide  consistently 
good  results  for  these  two  quantities.  Hierefore,  only  a,  e,  1,  and  a 
will  be  considered  In  presenting  the  reaialnlng  DCPLOC  reductions.  Hie 
observations  recorded  for  the  Fort  Sill,  Forrest  City  cooplex  are  plotte- 
In  Figure  12  for  six  revolutions  cf  Discoverer  XI  Including  172,  the  las., 
known  revolution  of  this  satellite.  Hie  DQPLOC  detemlned  position  for 
this  pass  Indicated  an  altitude  of  82  miles  as  '<die  satellite  crossed  the 
base  line  55  miles  west  of  Forrest  City.  A  cosiparlson  of  the  Space  Track 
results  with  the  DQPLOC  reductions  for  these  observations  Is  presented 
In  Figures  13  through  l6.  In  addition,  DQPLOC  reductions  have  been 
Included  for  three  revolutions  In  which  the  receivers  at  Forrest  City  and 
Aberdeen  Proving  Qround  tracaed  the  alr-bome  transmitter  In  the  satellite. 
In  Figures  17  through  20  a  similar  conparlson  Is  presented  for  six  separate 
passes  of  Transit  IB.  In  these,  all  observations  consist  of  data  obtained 
by  receivers  at  Forrest  City  and  Aberdeen  Proving  Ground  while  tracking 
the  on-board  transmitter. 

Finally,  In  Figure  21,  results  are  tabulated  for  a  reduction  based 
on  only  seven  frequency  observations.  Hiese  have  been  extracted  from  the 
cooqplete  set  of  observations  previously  presented  for  Sputnik  HI.  Hiey 
were  selected  to  serve  as  a  crude  example  of  the  type  of  reduction  required 
for  the  proposed  DQPLOC  scanning-beam  system.  Hie  exaisple  shows  that  the 
method  Is  quite  feasible  for  use  with  periodic,  discrete  measurements  of 
frequency.  Of  course,  the  proposed  system  would  normally  yield  several 
more  observations  than  were  available  in  the  exasQile. 

It  Is  noteworthy  that  numerous  solutions  have  been  obtained  with 
field  data  from  a  single  receiver  during  a  single  pass  of  a  satellite. 
Further,  these  measurements  have  been  confined  to  three  short  periods  of 
observatlcxi  within  a  two  to  three  minute  Interval.  Additional  receivers 
spread  over  greater  distances  would,  of  course,  considerably  enhance  the 
accuracy  of  the  results.  For  example,  a  system  with  two  receivers  and  a 
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IONOSPHERIC  MEASUREWn  rS 


Oue  of  the  methods  used  for  atudyleg  the  lonos^ere  Is  derived 
frost  the  mensureaent  of  Its  effect  oa  radio  waves  propagated  through  It. 

One  effect  observed  Is  the  Increase  in  propagation  phase  velocity  above 
that  occurring  In  free  space  (i.e.  the  wavelength  Is  longer  In  the  medlUB 
than  In  free  space).  At  the  frequencies  used,  30  to  300  me.,  the  index 
of  refraction  can  be  expressed  as 

40.25  K 

u  4  i - ,  (13) 

f 

o 

where  N  is  the  electron  density  In  electrons  per  cubic  aster  and  f  Is 
e  o 

the  transalttt^d  frequency.  Thus,  the  Index  of  refraction  Is  less  than 
unity  In  an  Ionized  aedlia  and  the  amount  by  which  the  Index  of  refraction 
Is  changed  is  Inversely  proportional  to  the  square  of  tbs  frequency. 

Use  Is  made  of  the  dependence  of  the  refractive  Index  ca  frequency  to 
detenlne  Ionosphere  electron  content.  Data  with  which  to  aaks  tha  electron 
content  computatlcn  are  obtained  by  Measuring  the  Doppler  frequency  shift 
on  two  baraonlcally  related  signal  frequencies  trsnsalttsd  from  s  satellite. 
The  Doppler  frequency,  sdperlapossd  upco  tbs  lower  of  the  two  frequencies, 
is  aultlpUed  by  the  ratio  of  tbs  two  signal  frequanclss  Miployed;  and  then 
the  result  la  subtracted  from  tbs  Doppler  observations  of  ths  higher  fre¬ 
quency  signal  to  yield  dispersive  Doppler  data.  The  dispersive  Doppler 
frequency  Is  proportional  to  the  time  rate  of  change  of  the  total  clecti"r* 
content  along  the  propagation  path, 

'dd  “  n  /V  '  (“> 

where  J  N  dr  is  the  total  electron  content  along  the  prcq^agallon  path,  r; 
o  * 

mus  a  measure  of  th«  cliange  in  electron  content,  over  a  time  Interval  of 
Interest,  Is  obtained  by  Integration  of  the  dispersive  Doppler  frequency . 
which  la  slJiply  a  cyc?c  coxmt  over  the  specified  interval. 
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Instrumentation  has  been  develcjied  for  measuring  and  recording 
dlsperclve  Doppler  frequency.  Satellites  carrying  radio  transmitters 
whose  frequencies  are  harmonically  related  serve  as  signal  sources.  In 
order  to  receive  signals  frca  satellites  at  great  distances  and  provide 
output  data  of  high  quality  It  Is  necessary  to  extremely  sensitive 
receiving  systems.  Narrow  bandwidth,  phase-locked,  tracking  filters  are 
used  to  provide  essentially  noise-free  Doppler  frequency  data.  Two 
channels  are  used  to  receive  the  two  'laimonlcally  related  signals.  The 
multiplication  of  the  Doppler  frequency  on  the  lower  frequency  signal  takes 
place  at  the  output  of  the  tracking  filter  so  as  not  to  degrade  the  signal- 
to-nolse  ratio.  Frequency  multiplication  prior  to  the  final  narrow-band- 
wldtii  filter  would  seriously  degrade  the  slgnal-to-nolse  ratio  of  the 
Doppler  signal.  A  special  broeuiband  frequency  multiplier  (ref.  4)  has  been 
developed  for  multiplying  audio  freqtiencles .  Ibe  technique  developed  Is 
unique  In  that  It  achieves  multiplication  of  an  audio  frequency  Doppler 
signal,  which  varies  many  octaves,  but  maintains  a  sinusoidal  output  wave¬ 
form.  The  multiplication  factor  can  be  any  product  of  two’s  and  three's 
(l.e.  2,  3,  4,  6,  8,  9,  — .).  This  frequency  multiplier  Is  basically  a 
combination  of  an  aperiodic  frequency  doubler,  push-push,  circuit  and  a 
bridge  configuration  trlpler  circuit.  Auxiliary  circuits  with  functions 
of  automatic  gain  control;  clipping,  differentiation  and  phase-locked 
tz*acklng  filtering  make  possible  a  sinusoidal  output  waveform. 

Figure  22  shows  a  block  dlagiam  of  a  receiving  system  for  Ionospheric 
measiurements  using  the  broeuiband  frequency  multiplier  (ref.  3).  Dispersive 
Doppler,  Faraday  rotation  and  satellite  ixstatloa  effects  on  the  signal  can 
be  separated  automatically  and  directly  recorded  as  r.hown  In  Figure  23. 

•nils  Is  a  portion  of  a  record  from  an  upper  atmosphere  sounding  rocket 
flight  In  which  a  two  frequency  trsnsmlcter  was  carried. 

Dispersive  Doppler  data  recorded  In  a  form  similar  to  that  shown  In 
Figure  23  can  be  counted  to  an  accuracy  of  -  0.1  cycle.  The  total  electron 
content  can  be  expressed  In  terms  of  dispersive  D->ppler  cycles  as 


(15) 


M^dr  =  (0^-K55j^) 


(K^  -  1)  15. X  10“ 


where  (0^  -  K0^)  is  the  integrated  dispersive  Doppler  frequency,  1?^  is 
the  lower  transmission  frequency,  Pg  is  the  higher,  K  is  the  ratio 
Consider  the  Transit  satellite  (i960  Eta),  where  »  54  me  and  P^  =■  324  me. 
The  Incremental  change  in  total  electron  content  for  each  dispersive 
Doppler  cycle  is 


3.24  X  10 


>  6.9  X  10^5  electrons 

square  meter  '  ' 


Therefore,  ti.e  covmtlng  accuracy  of  -  0.1  cycle  represents  a  meaaureing 

12 

sensitivity  to  the  change  in  ionosphere  electron  content  of  6.9  x  10 
electrona/square  meter.  This  sensttlvlty  is  high  enough  to  detect  small 
irregularities  in  the  ionosphere.  A  plot  of  dispersive  Doppler  data  and 
Integrated  dispersive  Doppler  frequency  is  shown  in  Figure  24  for  a  pass 
of  the  Transit  satellite  (i960  Eta)  on  17  November  I96O.  Irregular  hori¬ 
zontal  gradients  in  the  ionosphere  are  clearly  shown  by  the  variations  in 
the  dispersive  Doppler  frequency  that  are  evident  during  the  second  half 
of  the  satellite  pass.  This  cxirre  normally  has  a  relatively  smooth  "S" 
shape  under  undisturbed  geomagnetic  conditions.  It  is  of  considerable 
Interest  to  note  that  this  record  was  made  following  the  period  of  an 
extremely  severe  geomagnetic  disturbance.  Severely  disturbed  radio  cou- 
dltlons  existed  from  November  12  through  I8.  One  of  the  most  active  solar 
regions  observed  in  recent  years  was  reported  by  the  North  Atlantic  Radio 
Warning  Service  of  the  National  Bureau  of  Standards.  The  A- index  (a 
measure  of  geomagnetic  activity)  on  November  13  was  280,  the  hi  ..heat  recorded 
in  this  solar  cycle.  An  A- index  of  25  is  considered  a  disturbed  condition, 
therefore  260  represents  an  extremely  disturbed  condition.  An  unusuali.y 
high  magnetic  field  intensity  was  recorded  at  the  jelllstic  Research 
Laboratories  magr.ei'O'^c'cer  station  on  November  12  and  13  which  is  •'  m  ''n 
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Figure  Biufl  the  qu&lltatlve  agreeaent  between  the  irreguiariy  shaped 
dispersive  Doppler  curve  In  Figure  2k  and  the  disturbed  Icoosphere  Is  well 
established. 

?^e  Fhreulay  rotation  effect  can  also  be  used  to  detentlne  total 
Ionosphere  content.  Techniques  have  been  developed  for  separating  Faraday 
rotation  effects  from  satellite  rotation  effects  by  the  uae  of  opposing 
circularly  polarized  antennau  and  a  sequence  of  electronic  alxera  aa  ahcnm 
In  Figure  22.  When  the  satellite  spins  very  slowly,  a  simpler  method  of 
determining  Faraday  rotation  cycles  by  counting  received  signal  aa^pUtude 
nulla  can  be  used.  A  linearly  polarized  receiving  antenna  is  used  In  this 
case.  A  plot  of  Ionospheric  electron  content  Is  shown  In  Figure  26,  obtained 
by  losing  received  signal  amplitude  null  data  from  a  pass  at  the  Transit 
satellite  (i960  Bta).  'Rje  computation  methods  of  BowhlU  (ref.  6)  and 
Garnott  (ref.  7)  were  used  In  the  earliest  studies.  A  coaqdete  ray 
tracing  program  baaed  on  that  of  Little  and  Lawrence  (ref.  8)  Is  In  prepa¬ 
ration  to  provide  nore  accuracy  and  to  eliminate  several  assumptions  and 
restrictions  of  the  ecu:ly  methods. 
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FIGUHE  26 


CQHCLUSIQNS 


The  Infonaatlon  on  the  Ionosphere,  o'  sained  hy  the  methods  des- 
erlhcd,  makes  It  possible  to  correct  refra'ition  errors  and  obtain  more 
accurate  orbital  parameters  frot  Doppler  lata.  An  interesting  exas^ple  of 
the  Ionospheric  effect  on  orbital  accuracy  was  ooserred  in  the  coQputatlon 
of  the  orbital  parameters  shown  in  Figures  17>  l8,  19,  and  20.  The 
computation  was  first  attempted  using  the  complete  "S”  curve  including 
the  relatively  constant  frequency  limbs.  Bie  limbs  represent  data  ob¬ 
tained  during  the  emergence  of  the  satelT ite  from  the  horizon  and  recession 
Into  the  horizon.  The  orbit  obtained  was  appreciably  different  from  that 
published  by  Space  Track.  Another  eoaputatlon  was  made  using  only  the 
center  portion  of  the  "S"  curve,  while  disregarding  the  limbs.  The  so¬ 
lution  was  sireatly  improved  and  the  results  agreed  very  well  with  Space 
Track  data.  This  points  out  the  large  refractive  effect  the  ionospbere 
Introduces  at  low  elevation  angles  of  transmlaalon.  Fortunately,  an 
orbital  solution  can  be  computec.  from  Doppler  data  obtained  at  quite  high 
elevation  angles,  thereby  minimizing  the  refractive  error-. 

A  program  haa  been  initiated  to  ccablne  D(^ppler  frequency  observations 
with  electron  content  data  in  an  i-terative  coerputlng  process  to  Increase 
the  accuracy  of  toe  orbital  deteiainatlon.  the  computation  will  be  initi¬ 
ated  by  detenlnlng  an  orbit  from  the  uncorrected  Dopplar  observations. 

The  electron  content  data  and  this  approximate  orbit  will  be  combined  tc 
compute  corrections  for  the  original  Doppler  frequency  meaaurements . 

Using  -the  latter,  the  process  will  be  iterated  until  the  refractive  error 
hac  been  minimized  in  tbe  Doppler  data  and  hence,  in  the  coaptited  orbital 
parameters  as  well. 


IF.B. 

R.  B.  PATTOf,  JR*' 

V.  W.  RICuARD 
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to  Kioialu  nodoB  error*  la  oiftler  freiji»?=oi-  .^sOaraKiit*  derieed  fro*  ve«>  to  BlnlKlte  rudoa  error*  la  Doppler  frequency  aeatureBente  derived  froB  weak 

•l^iale  trBiaalttad  over  eetreae  ntn«ra.  alsnale  traniBltted  over  eztreBC  raii«ea. 

tie*  la  alao  Bade  of  aac.  Uitf  .  -ppier  da*  to  deteiBlnc  t>  •  total  clectroa  Uae  la  alao  aade  of  aatelllte  Doppler  daU  to  deteiBlne  the  total  electron 

coateo*.  In  the  toBoaptaerc-  a  tect-'  .ue  for  aoazurlac  Uie  dlep  -alee  Don>l*r  content  la  the  loooaphere.  A  technique  for  Beaaurtng  the  dlcperalve  Doppler 

effect  li  deaerthad.  fleaul.ta  of  coapj'aUon*  with  actual  fl.-lu  data  are  preaeoted.  effret  la  deacrlbed.  Beaulta  of  coBputatlona  with  actial  field  data  are  preaented 


